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We consider theoretically as a function of temperature the plasmon mode arising in three- 
dimensional Dirac liquids, i.e., systems with linear chiral relativistic single-particle dispersion, within 
the random phase approximation. We find that whereas no plasmon mode exists in the intrinsic 
(undoped) system at zero temperature, there is a well-defined finite-temperature plasmon with 
superlinear temperature dependence, rendering the plasmon dispersion widely tunable with temper¬ 
ature. The plasmon dispersion contains a logarithmic correction due to the ultraviolet-logarithmic 
renormalization of the electron charge, manifesting a fundamental many-body interaction effect as 
in quantum electrodynamics. The plasmon dispersion of the extrinsic (doped) system displays a 
minimum at finite temperature before it crosses over to the superlinear intrinsic behavior at higher 
temperature, implying that the high-temperature plasmon is a universal feature of Dirac liquids 
irrespective of doping. This striking characteristic temperature dependence of intrinsic Dirac plas- 
mons along with the logarithmic renormalization is a unique manifestation of the three-dimensional 
relativistic Dirac nature of quasiparticle excitations and serves as an experimentally observable 
signature of three-dimensional Dirac materials. 
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In Dirac semimetals the valence and conduction bands 
touch only in isolated points of the Brillouin zone with a 
Fermi energy that is tuned close to these band touching 
points, giving rise to particle excitations with an effec¬ 
tive linear relativistic band dispersion £±(k) = ±hvpk. 
Much of the interest in these novel materials stems from 
the topological properties of Weyl semimetals which are 
manifest in anomalous transport properties mm and 
Fermi arc surface states |3 . Recent hallmark experi¬ 
ments have realized these materials in CdsAs 2 , NasBi, 
and TaAs and have demonstrated the Dirac cone struc¬ 
ture in ARPES measurements Hi] and detected topo¬ 
logical surface states in NasBi and TaAs |9j[l0]. These 
initial experiments, which serve to verify the expected 
band structure, motivate the search for definitive ex¬ 
perimental signatures of Dirac materials which establish 
unique observable effects distinct from standard electron 
gases and normal metals. In this Rapid Communica¬ 
tion, we show that the collective charge oscillation of a 
Dirac liquid, the plasmon, is widely tunable as a super- 
linear function of temperature (with ^ —T/lnT), which 
allows us to extract high-precision information on the 
Dirac material. Moreover, the dispersion reveals a subtle 
logarithmic electron charge renormalization effect which 
is an additional specific signature of a three-dimensional 
Dirac system. Note that previous work exclusively stud¬ 
ies the plasmon mode at zero temperature mHE]. In this 
Rapid Communication, we establish a different mecha¬ 
nism for finite-temperature plasmons which arises due to 
the thermal excitation of electrons from the valence to 
the conduction band. It turns out that this interband 
mechanism overwhelms the low-temperature canonical 
intraband mechanism at sufficiently large temperature, 
resulting in a universal behavior that is independent of 


initial detuning. 

The collective plasmon mode of an electron system is 
defined as the complex pole z = uj(q) — ry(g) of the di¬ 
electric function e(z, q); i.e., it solves 

s(z, q) = 1 - Vqll(z, q) = 0, (1) 

where n (z, q) denotes the polarizability and Eq is the 
Coulomb interaction E q = 47re 2 /ftq 2 with e being the 
electron charge and k, the effective dielectric constant of 
the medium. The dimensionless interaction strength (the 
Dirac fine structure constant) is defined by the ratio of 
interaction and kinetic energy 


Here, because of linear band dispersion, we have the spe¬ 
cial situation that both kinetic and potential energy scale 
as n 1 / 3 with carrier density, implying that a is density 
independent. This manifests a scale invariance of the the¬ 
ory on the classical level. By contrast, the nonrelativistic 
system with quadratic dispersion h 2 k 2 /2m is generically 
scale-invariant only for 1/r 2 long-range interactions (or, 
for short-range interactions, by fine-tuning to the unitary 
limit of infinite scattering length), implying that for the 
nonrelativistic Coulomb interactions the effective inter¬ 
action strength r s ~ n -2 / 3 can be changed by altering 
the density. 

The form of the plasmon dispersion both at zero and 
finite temperature is determined by dimensional analy¬ 
sis: at zero temperature and finite density, the plasmon 
dispersion uj p ~ e F ^ n 1 / 3 , where ep is the Fermi energy 
and n is the doping-induced carrier density, approaches 
a constant as is characteristic for three-dimensional sys¬ 
tems. It should be noted that in contrast to the electron 
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liquid the plasmon frequency uo p of the Dirac liquid is an 
entirely quantum mechanical quantity m with h appear¬ 
ing explicitly in the leading-order low-momentum term of 
the plasmon dispersion, uo p ^ 1 /\/h [whereas the electron 
gas with uj p = A.'Kne 2 jKm receives quantum corrections 
only in higher order 0(q 2 /h 2 )\. At zero temperature and 
zero density (the intrinsic case), by dimensional analysis 
we expect the plasmon dispersion to vanish. At finite 
temperature, on the other hand, thermal excitations of 
particles from the valence to the conduction band are 
not gapped and a solution of Eq. © exists even at zero 
doping. Here, the low momentum plasmon frequency 
should be proportional to temperature, uj p ~ T. It turns 
out that this simple scaling analysis is not quite correct 
and the full solution predicts uj p ~ —T/lnT: both plas¬ 
mon dispersion and damping receive a quantum correc¬ 
tion due to the renormalization of the Coulomb interac¬ 
tion, a field-theoretical effect in interacting many-body 
systems which does not possess an analog in standard 
electron gases with parabolic band dispersion or Dirac 
materials in other dimensions m- 

In this Rapid Communication, we formalize this dis¬ 
cussion by a full random phase approximation (RPA) 
calculation of the plasmon dispersion at finite temper¬ 
ature and density. We compute analytically the plasmon 
dispersion and damping, finding a characteristic superlin- 
ear temperature dependence which includes a logarithmic 
renormalization correction for the intrinsic Dirac liquid. 
For the extrinsic doped system, our calculations indicate 
a crossover from a low-temperature extrinsic regime to a 
high-temperature intrinsic regime where the plasmon be¬ 
haves as in the undoped system. We also complement all 
our analytical findings with full numerical calculations 
of the plasmon properties. In addition to the plasmon 
dispersion, we present numerical results for the dielec¬ 
tric function as a function of momentum, frequency, and 
temperature for the Dirac-Weyl system. 

Within RPA, we compute the plasmon mode in Eq. 0 
using the noninteracting polarizability 


n (*><o = ££ 

kss' 


/( gj (k))-/( gj ,(kQ) 

hz + £ s ( k) — S s > (k 7 ) 


F SS '( k,k'). 


( 3 ) 


Here, k' = k + q, f(s s ( k)) = np(£ s (k) — g) is the 
Fermi-Dirac distribution, g the chemical potential, and 
F ss /(k, k') = (1 + ss' cos#)/2 is the overlap of two wave 
functions with energy e s (k ) and e s f ( k ') which can be eval¬ 
uated in terms of the relative angle 0 between k and kb g 
accounts for a degeneracy of Dirac cones in the material. 
Equation © diverges logarithmically in momentum na- 
IDj, and we introduce a scale A beyond which we cut off 
the linear Dirac dispersion. Throughout this Rapid Com¬ 
munication, we set ks = 1 but retain the dependence on 
h where appropriate. 

We can split the polarizability into a zero-temperature 
intrinsic contribution which encompasses the complete 


cutoff dependence and a correction due to finite tem¬ 
perature and doping. At zero temperature, the intrinsic 
polarizability reads 


n(w,q) = - 


gq 


2AlT 2 hVF 
2 

Q(uj 
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uS 2 /v 2 f — q 2 
v F q). 


( 4 ) 


Substituting this in Eq. 0 , we notice that the expression 
does not involve a and A separately but combined in the 
form of a single dimensionful scale A l = Ae 3?r /^ a , the 
Landau pole m ■ In typical Dirac materials with k ^ 1, 
we find the following hierarchy of scales: 


kp <C A <C A^. 


( 5 ) 


Note that even a weak screening increases the Landau 
pole scale exponentially compared to the cutoff A. This 
appearance of a dimensionful scale is known as dimen¬ 
sional transmutation and signifies the breakdown of clas¬ 
sical scale invariance. This so called scale anomaly, 
where a system that is scale-invariant on a classical 
level acquires a dimensionful scale through the process 
of renormalization, is present in other condensed mat¬ 
ter or atomic systems as well: the most pivotal exam¬ 
ple is the standard BCS theory with attractive contact 
interaction of strength g. Here, the renormalization of 
the dimensionless coupling implies that the combination 
A = c or>e l / VN ^ 9 being the Debye frequency, i.e., 

the cutoff scale of the theory) is an observable quantity 
independent of the renormalization group scale. More¬ 
over, in two-dimensional atomic Fermi gases, the break¬ 
ing of scale invariance is conjectured to cause a shift 
in the collective breathing mode 2D 122] and to induce 
logarithmic scaling violations in the high-frequency tail 
of radiofrequency absorption spectra [23] and other re¬ 
sponse functions [24]. It should be noted, however, that 
in the above examples the scale arising due to dimen¬ 
sional transmutation is an inherently low-energy scale 
(the superconducting gap or the dimer bound state en¬ 
ergy, respectively), whereas in the Dirac materials dis¬ 
cussed in this work the Landau scale is larger than the 
cutoff A. We remark that this scale anomaly is inherently 
a nonperturbative effect involving an essential singular¬ 
ity in the coupling constant a. Note that the cutoff A 
corresponds to the momentum scale beyond which the 
dispersion is no longer linear. For Weyl materials with 
Weyl pairs in close proximity in the Brillouin zone (such 
as the pyrochlore iridates or TaAs mm), the Landau 
pole could be sufficiently small to ensure large logarith¬ 
mic scaling corrections. 

In the following, we analyze the plasmon dispersion 
and dielectric function at finite temperature both ana¬ 
lytically and numerically. We consider first the intrinsic 
and then the extrinsic plasmon. For intrinsic Dirac ma¬ 
terials, the polarization at small momentum q and high 
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FIG. 1. Intrinsic finite temperature dielectric function with 
KvfAl/T — 5 • 10 3 (A /kF — 100 and ga — 2.4). The scale 
is logarithmic. The white dashed line indicates the plasmon 
dispersion. 


temperature T reads 
nr ) = gq2 i 

' ^ 127 t 2 Hvf 7t e~ 1E T 18H 3 vf w 2 96ttvfT' 

(6) 

where 7 # is the Euler constant, giving rise to the plasmon 
dispersion 


fuj p = 7 tT a 


31n/mpA£/7re 


l27rgaT 
k(T) 3' U 


We have expressed the frequency in terms of the Landau 
scale Al as well as the more conventional form with an 
effective temperature dependent dielectric constant 


^= 1+9 £ ln 


hvp A 

7r e~^ E T 


ga , HvfAl 

— la-— • 

3tt Tre 7^T 


(8) 


The plasmon temperature dependence is of superlinear 
form l d p ~ — T/lnT, which is one of the main results 
of this Rapid Communication. For T = 100 K and 
HvfAl = 500 K and 5000 K, for example, the plasmon 
energy 0 is fujp = 20 meV and lOmeV, corresponding 
to a THz frequency range. Note that there is a significant 
depression of the plasmon frequency compared to a stan¬ 
dard Dirac plasma [i.e., a material without a hole band 
of negative energy], in which renormalization effects are 
absent and k(T) = 1. Both the anomalous temperature 
dependence and the change in the plasmon frequency we 
predict here should be readily observable. 

Expanding the polarizability to leading order around 
we obtain the damping of the plasmon mode 


fry P = 


_7r^T_ 

24 In 2 hv F h. L /ire-7ET 


( ga \ 2 7 tT 

\6k(T) J 6 


(9) 




FIG. 2. Extrinsic 3D plasmon dispersion and damping for 
ga — 2.4 and A/fcp = 100 (corresponding to A l/^f — 5 • 
10 3 ) for various momenta (bottom to top) q/kF — 0,0.3, 0.6, 
and 0.9. The red continuous lines indicate the exact results 
derived in Eqs. |7|, ([ 9 ]), |l0| , and (14). 


The plasmon mode is only well defined for 7 p /cj p <C 1. 
While both damping and dispersion increase with tem¬ 
perature, it turns out that the ratio 

ip_ _ _^_ = jl ay — r 

^p HvfAf/ 7re~^ E T 72 y 2 ^^(T 1 )^ 

( 10 ) 

is numerically small. Indeed, the renormalization of the 
electron charge with k(T) > 1 decreases this ratio and 
strongly suppresses the damping relative to the frequency 
for T < HvfAl- This is another central result of this 
work. To further illustrate this point, we show in Fig. [l] 
the full numerically computed intrinsic dielectric func¬ 
tion for a fixed Landau pole scale HvfAl/T = 5-10 3 . 
It is apparent that density fluctuations are dominated 
by a sharp plasmon dispersion as discussed above. The 
plasmon mode at small momentum starts off at constant 
frequency given by Eq. (|7|) with a parabolic dispersion, 
which exists up to large momenta and eventually merges 
with the particle-hole continuum. 

The extrinsic plasmon displays a crossover from a gen¬ 
uine extrinsic behavior governed by thermal intraband 
excitations to the intrinsic behavior at high temperature. 
In contrast to the intrinsic case, the chemical potential 
here is nonzero and depends on temperature. It is fixed 
by the relation 

|pE n +( k ) + _ !] = n ’ ( n ) 

k k 

where n s are the Fermi-Dirac distributions with energy 
£ s (k) = shvFk and n = gk^/Gi: 2 is the zero temperature 
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carrier density. We find that the chemical potential has a 
quadratic deviation from sp at low temperature fije F — 
1 — 7r 2 T 2 /3T| 1 , where Tp denotes the Fermi temperature, 
and goes to zero as a power law at high temperature 
g/e F = T 2 /n 2 T 2 . 

At zero temperature, the extrinsic plasmon frequency 
is given by HHIE1IEI 


^ = \/i^wS7 


£p = 


O7T/^0 


( 12 ) 


where is a renormalized dielectric constant 

ga A l A 

fto = — In —— = 1 H-In ——. 

37r 2 kp 37r 2 kp 


(13) 


Note that the extrinsic plasmon frequency goes as the 
cube root of the doping density up to logarithmic correc¬ 
tions. There are two low-temperature corrections to the 
plasmon dispersion: first, the temperature dependence of 
the chemical potential discussed in the previous section, 
and second, a Sommerfeld correction due to the smooth¬ 
ing of the Fermi surface. This gives 


co p (T) = 




0(T 4 /T£). 


(14) 


The damping of the plasmon mode at low temperatures 
requires the thermal excitations of particle-hole pairs and 
is therefore suppressed by powers of . 

The intrinsic behavior discussed above also describes 
the extrinsic material at higher temperatures 0(Tp ) as 
can be seen by computing the high-temperature limit of 
the extrinsic polarizability, 


n(cj,g) = - 


gq 


In 


127 X 2 hvp 

1 - 


igq 2 uj 


967 xvpT 


IwfA gifr 2 

ne~^ E T 6Tr 2 h 3 VFU 2 
H 2 lo 2 g 2 
48T 2 ~~ 4T 2 ’ 


7T 2 g 2 

-I- — 

3 T 2 


(15) 


implying that corrections to the intrinsic plasmon dis¬ 
persion at high temperature are of order 0(g 2 /T 2 ) = 
0(1/T 6 ) and thus highly suppressed. 

Equations (14) and (15) imply that the plasmon fre¬ 
quency displays a nonmonotonic behavior where it first 
decreases with increasing temperature manifesting the 
extrinsic behavior, reaches its minimum, and then in¬ 
creases with the characteristic intrinsic superlinear be¬ 
havior. The minimum, which is reached at a temper¬ 
ature scale O(Tp), marks the crossover between gen¬ 
uine extrinsic and intrinsic behavior. The measurement 
of the plasmon dispersion as a function of temperature 
thus provides an independent way to determine the ex¬ 
trinsic Fermi energy. The analytic results for the plas¬ 
mon crossover derived in the previous section are in close 
agreement with a full numerical calculation of the plas¬ 
mon dispersion as shown in Fig. [2] The analytical results 


in Eqs. 0,0, ©, and ( p~4| ) are indicated by red lines. 
It is apparent that the full result interpolates closely be¬ 
tween the low-temperature and the high-temperature ex¬ 
act results. Here, we have the interesting situation of 
a rapid crossover between the low-temperature regime, 
where density fluctuations are completely determined by 
the intraband excitations of the extrinsic system [up to 
a renormalization of the effective coupling /c(fcp)], and 
the high-temperature regime completely characterized by 
interband excitations. It is important to note that ir¬ 
respective of the initial detuning, any extrinsic effects 
will be washed out at sufficiently high temperature. Fig¬ 
ure [2] also shows the plasmon frequency and damping 
at nonzero momentum where no analytical results ex¬ 
ist. While the damping is almost unchanged in the mo¬ 
mentum range considered here, the plasmon dispersion 
increases with increasing momentum. At high temper¬ 
ature, the intrinsic plasmon dispersion persists (with a 
higher offset compared to the zero-momentum case), but 
the minimum is less pronounced and eventually vanishes. 

In the following, we discuss the crossover from extrinsic 
to intrinsic behavior for the full dielectric function. Con¬ 
venient for our numerical computations is the following 
integral representation for n(cj, q) at finite temperature: 


n(w,q,/x, 



dg! 


n(u,q, g',T = 0) 

4Tcosh 2 ((/z' — g)/2T) ’ 


(16) 


where the full extrinsic polarizability can be evalu¬ 
ated in closed analytical form at zero temperature HU. 
The crossover from extrinsic low-temperature to intrin¬ 
sic high-temperature behavior is apparent in the full di¬ 
electric function as can be seen in Fig. [3] While at zero 
temperature, interband excitations only exist for frequen¬ 
cies uj > 2kp — q, this sharp boundary is quickly washed 
out even for moderate temperatures. At higher tempera¬ 
tures, the dielectric function is equal to the fully intrinsic 
result and corresponds to Fig. 0 up to a rescaling of 
units. 

The discussion so far has focussed on the case of com¬ 
pletely symmetric Dirac cones, i.e., with single-particle 
dispersion e s (k) = shvpk. By contrast, all experi¬ 
mentally realized Dirac materials (and indeed all the¬ 
oretically conjectured systems), have a strongly asym¬ 
metric structure which is reflected in the particle dis¬ 
persion e s (k) = shyjvF X k% + Vpyky + v 2 Fz kl with typ¬ 
ically vp z vpx^vpy. We note that our discus¬ 
sion fully applies to these systems by a simple replace¬ 
ment of the Fermi velocity v F by the geometric mean 
[■ vpxVpyVpz ] 1//3 and rescaling the momentum q jvp 
{qx/vFx,qy/vFy,q z /vFz) in all of our results. We also 
note that although our results are obtained within RPA, 
they should be essentially exact for realistic Dirac-Weyl 
systems where the typical degeneracy factor g 1 by 
virtue of the many-valley multiple Dirac cone nature of 
these systems [g = 24 in the recently discovered Weyl- 
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FIG. 3. Extrinsic 3D finite temperature dielectric function for ga — 0.62 and A/fcp = 100 and various temperatures T/Tf = 
0.01,0.05,0.1,0.15,1., and 2. The white dashed line indicates the plasmon dispersion. The crossover from the extrinsic zero- 
temperature dispersion to the intrinsic high-temperature limit is clearly visible. 


semimetal TaAs wm and all higher-order many-body 
corrections are suppressed by higher powers of 1/g. The 
logarithmic ultraviolet renormalization is unaffected by 
higher-order corrections by virtue of the renormalizabil- 
ity of the Coulomb interacting Dirac systems. 

We have presented a comprehensive study of collec¬ 
tive plasmon modes in 3D Dirac materials at finite den¬ 
sity and temperature. We find that the plasmon mode 
is extremely well defined with damping rates more than 
an order of magnitude smaller compared to the energy. 
Moreover, the plasmon dispersion is super linear in tem¬ 
perature with logarithmic Landau pole corrections due 
to the renormalization of the electron charge. Our re¬ 
sults are of fundamental interest as a manifestation of 
charge renormalization similar to the one found in quan¬ 
tum electrodynamics, and of significant practical rele¬ 
vance as temperature offers a way to tune the plasmon 
dispersion at will. Our predictions should be directly ob¬ 
servable in present-day experiments and should provide 
definitive signatures for the Dirac nature of candidate 
electron systems. 

Note added. Recently, Ref. [26] appeared, which 
studies intrinsic finite-temperature plasmons in Dirac 
semimetals and reports the intrinsic plasmon dispersion. 
Our result reported in Eq. 0 differs from Ref. m by a 
factor of 7re _7E and a different definition of the Coulomb 
interaction, which includes a dielectric constant k. 

This work is supported by LPS-MPO-CMTC. 
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